Mixtures of ultracold fermions with unequal masses 
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We analyze the phase diagram of superfluidity for two-species fermiou mixtures from the Bardeen- 
Cooper-Schrieffer (BCS) to Bose-Einstein condensation (BEC) limit as a function of scattering 
parameter, population imbalance and mass anisotropy. We identify regions corresponding to normal, 
or uniform/non- uniform superfluid phases, and discuss topological quantum phase transitions in 
the BCS, unitarity and BEC limits. We derive the Ginzburg-Landau equation near the critical 
temperature, and show that it describes a dilute mixture of paired and unpaired fermions in the 
BEC limit. We also obtain the zero temperature low frequency and long wavelength collective 
excitation spectrum, and recover the Bogoliubov relation for weakly interacting dilute bosons in the 
BEC limit. Lastly, we discuss the effects of harmonic traps and the resulting density profiles in the 
BEC regime. 
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I. INTRODUCTION 

The evolution from Bardeen-Cooper-Schrieffer (BCS) 
state to Bose-Einstein condensation (BEC) is an im- 
portant topic of current research for condensed mat- 
ter, nuclear, atomic and molecular physics communi- 
ties. Recent advances in atomic physics have allowed 
for the study of superfluid properties in symmetric two- 
component fermion superfluids (equal mass and equal 
populations) as a function of scattering length, where 
the theoretically predicted crossover from BCS to BEC 
was observed Jj, Q ■ 

Since the population of each component as well as the 
interaction strength between two components are exper- 
imentally tunable, these knobs enabled the study of the 
BCS to BEC evolution in asymmetric two-component 
fermion sijperfluids (equal mass but unequal popula- 
tions) 0, Si- In contrast with the crossover physics 
found in the symmetric case, these experiments have 
demonstrated the existence of phase transitions between 
normal and superfluid phases, as well as phase sep- 
aration between superfluid (paired) and normal (ex- 
cess) fermions as a function of pop ulation imbalance 

a S i, M, MMMM El, mM, MMMM^- 

Arguably one of the next frontiers of exploration 
in cold Fermi gases is the study of asymmetric two- 
component fermion superfluidity (unequal masses and 
equal or unequal populations) in two-species fermion 
mixtures from the BCS to the BEC limit. Earlier works 
on two-species fermion mixtures were limited to the BCS 
regime [1, 0, H, HI]. However, very recently, the evolu- 
tion from BCS to BEC was preliminarily addressed in ho- 
mogenous systems as a function of population imbalance 
and scattering length especially for ^Li and *'^K mixture 
and also for other mixtures including ^Li and ^^Sr and 
and ®^Sr mixtures [IS, HE HI. In addition, the su- 
perfluid phase diagram of trapped systems at unitarity 
was also analyzed as a function of population imbalance 
and mass anisotropy [13] ■ 

In this manuscript, we study the BCS to BEC evolu- 



tion of asymmetric two-component fermion superfluids as 
a function of population imbalance and mass anisotropy, 
and extend our previous works [24 . in this area. Our 
main results are as follows. For a homogeneous system, 
we analyze the zero temperature phase diagram for ultra 
cold fermion mixtures as a function of mass anisotropy 
and population imbalance. We identify regions corre- 
sponding to normal, uniform or non-uniform superfluid 
phases, and discuss topological quantum phase transi- 
tions in the BCS, unitarity and BEC limits. We de- 
rive the Ginzburg-Landau theory near the critical tem- 
perature, and show that it describes a dilute mixture of 
weakly interacting paired and unpaired fermions in the 
BEC limit. We also derive the zero temperature low fre- 
quency and long wavelength collective excitation spec- 
trum for zero population imbalance, and recover the Bo- 
goliubov relation for weakly interacting dilute bosons in 
the BEC limit. In addition, we describe analytically the 
phase separation boundaries of the resulting Bose- Fermi 
mixture of paired fermions and unpaired fermions in the 
BEC limit. Lastly, we discuss the effects of harmonic 
traps and the resulting density profiles of paired and un- 
paired fermions in the BEC regime. 

The rest of the paper is organized as follows. In Sec.lTTl 
we introduce the imaginary-time functional integration 
formalism, and obtain the self-consistency (order param- 
eter and number) equations. In Sec. Illli we discuss the 
evolution from BCS to BEC superfluidity at zero tem- 
perature within the saddle point approximation, and we 
analyze the order parameter, chemical potential, and sta- 
bility of the saddle-point solutions. In Sec. lIVi we discuss 
gaussian fluctuations near the critical temperature, and 
we obtain the low energy collective excitations and cor- 
rections to the saddle point phase diagram in the BEC 
region at zero temperature. In addition, we discuss the 
effects of harmonic traps on the density profile of paired 
and unpaired fermions at zero temperature. A summary 
of our conclusions is given in Sec.|Vl Lastly, we present in 
Appendices [XI [Bl and [Cl the elements of the inverse fluc- 
tuation matrix, and their low frequency and long wave- 
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length expansion coefRcients at finite and zero tempera- 
tures. 



II. FUNCTIONAL INTEGRAL FORMALISM 

To describe a dilute two-species Fermi gas in three di- 
mensions, we start from the Hamiltonian (?i = 1) 

H = Y^ ^k.^aL^flk,. + E '^(k, k')6^6k',q, (1) 



k.cr 



k,k',q 



where the pseudo-spin a labels both the type and hy- 
perfine states of atoms represented by the creation op- 
erator a^^, and b^^ = 4+q/2,T«-k+q/2,i- H*^^^' ^k,<T = 
ek,(T — /^cri where ek.o- = k"^ /{'2m„) is the energy and /i^- 
is the chemical potential of the fermions. Notice that we 
allow for the fermions to have different masses mg. and 
different populations controlled by independent chemi- 
cal potentials /i^- The attractive fermion-fermion inter- 
action ?7(k, k') can be written in a separable form as 
L/(k, k') = -gr*rk' where g > 0, and Tk = 1 for the 
zero-range s-wave contact interaction considered in this 
manuscript. 



Effective action 



The gaussian effective action for H is 



^Gauss = 5o-l-^^At(g)F-i(q)A(g), (2) 

where (3 — 1/T and q — (q, Wf) with bosonic Matsubara 
frequency = 2£Tr//3. Here, the saddle point action is 
given by 



So=P 



|Ao| 



^Trln[G^P(p)//3]-\ (3) 



where p = (k, Wi) with fermionic Matsubara frequency 
Wi = {2£ + l)7r//3, Ak — AoFk is the order parameter, 
and (G''P)~'^(p) is the inverse fermion propagator. The 
matrix {G^P)~^{p) is defined by 



(Ci^P)-^(n) - f *"^«~^k,T Ak 



(4) 



Furthermore, the vector A''{q) is the order parameter 
fiuctuation field and F~^{q) is the inverse fluctuation 
propagator. The matrix elements of F~-^{q) are given 

by 

Fm(9) = 7-^EGT^T(f +P)G!:'x(f "P)|rkP(5) 



(6) 



Notice that F^Hq) ^ (Frj)^?) and F^^^{q) = F^^^{-q). 
These matrix elements are described further in ap- 
pendix El The fiuctuation term in the action leads to 
a correction to the thermodynamic potential, which can 
be written as ricauss = J^o + ^^fiuct with fig = Sq / fi and 
rjfl^rt = (l//3)Eglndet[F~i((7)//3]. Next, we derive the 
self-consistency equations. 



B. Self-consistency equations 

The saddle point condition (55o/(5Aq = or the rela- 
tion 

Ak' = - ^ lim V C/(k, k')G!P, (p) exp(zw;,r) (7) 



leads to an equation for the order parameter 



2E] 



kH 



k, + , 



(8) 



where X\^^± = ('^k,i ± '^k,2)/2 with X^s = 
tanh(/3i<^k,s/2). Notice that, at low temperatures T « 0, 
0{—Ek,s) = lini/3^oo '^k.s, where 0{x) is the Heaviside 
function. Here, -Ek,± = (-E-k.i ± -E'k,2)/2 and ^k,± = 
(Ck,T ±^k,i)/2 = fcV(2'™±) -M±, where 

= (el+ + |Akn^/'+a,-, (9) 

i?k,2 = (ei+ + iAkni/2-a,_ (10) 

are the quasiparticle and negative of the quasihole ener- 
gies respectively, m± = 2m^mi/ {jni ± m|) and ^± = 
(^1 ± fJ.i)/2. Notice that m+ is twice the reduced mass 
of the "f and | fermions, and that the equal mass case 
corresponds to |m_| — > oo. As usual, we eliminate g in 
favor of the scattering length ap via the relation 



1 _ ni+V ^ y iFkp 



(11) 



where ek,± = (ek,T ± ek,x)/2. 

The order parameter equation has to be solved self- 
consistently with number equations — — 
which have two contributions iV^ = Nq ,, + -/Vfluct^a-- The 
first term Nq^^ = —dVlo/d^a or the relation 



0,cr 



r^O ^ — ^ 



WlT) 



(12) 



leads to the saddle point contribution, and is given by 

l-7cr'%'k,- Ck 



k 



2E, 



-X, 



k,+ 



(13) 



Here, 7| = +1 and 7^ = —1. Similarly, the second term 
A^fluct.o- = ^Srifiuct/c^Mo- leads to the fluctuation contri- 
bution, and is given by 



A"fluct,£T - -75E 



^[deiF-\q)]|^^l„ 



13^ detF-i(q) 



(14) 
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While the saddle point contribution is sufficient for a 
semi-quantitative analysis at zero temperature (T « 0), 
inclusion of the fluctuation contribution is necessary to 
recover BEC physics at finite temperatures (T Tc). 

When populations of the pseudo-spin components are 
balanced {N^ = Ni), the results for |Ao| and /Lt+ is 
irrelevant in this case) of m| ^ m| case can be obtained 
from the results for |Ao| and /i of m = m-f = case via 
a substitution of to — )■ to_|_. However, when populations 
of the pseudo-spin components are imbalanced {N^ ^ 
7V|), we need to solve all three self-consistency equations, 
since population imbalance is achieved when either -Bk.i 
or i?k,2 is negative in some regions of momentum space, 
which is discussed next. 



III. SADDLE POINT RESULTS 

At low temperatures T w 0, the saddle point self- 
consistency (order parameter and number) equations 
are sufhcient to describe the evolution of superfluid- 
ity from the BCS to the BEC limit. In this section, 
we analyze amplitude of the order parameter |Ao| and 
chemical potentials /io- as a function of mass anisotropy 
rrir = m^/mi and population imbalance P = iV_/iV+ 
for a set of fixed scattering parameters l/{kF,+aF)- Here 
N± = (iV| ± Ni)/2 and k^^^ = {k^,^ ± k%^^)/2. Be- 
cause of the parabolic dispersion relation, the density 
of t fermions is = A:|.|/(67r^) and the density of [ 
fermions is Ni — kp ^/ {6tt'^). Here, the Fermi momenta 
/cf.t ^-iid kF,i are determined from the Fermi energies 
eF,cT = (t/(2?t1(t). Therefore, the total fermion density 
is 

n = ni+ni = —^. (15) 

Using the notations described in the preceding para- 
graph, we can solve the self-consistency Eqs. ([8]), (fTT|) 
and For instance, in Fig.[Tl we show self-consistent 
solutions of |Ao|, /i+ and /i_ (in units of £f.+) as a func- 
tion of rrir when P — 0.5 and l/{kF,+aF) — 0. Here, 
ep^ = kp _^./ (2m±). With increasing mass anisotropy 
(or decreasing m^), we find that both |Ao| and ii+ in- 
crease until rrir ~ 0.4. However, further decrease in 
beyond to^ w 0.4 leads to a saturation of both |Ao| and 
with a small cusp in both quantities. The cusp is 
best seen in Fig. [T] at small grazing angles. Therefore, 
the evolution from rrir = 1 to to^ ^ is non-analytic 
when TOr « 0.4, and the evolution is not a crossover. 
These cusps in |Ao| and are more pronounced for 
higher |P|, and they signal a topological quantum phase 
transition discussed below. Notice that, for P — 0, the 
evolution of |Ao| and /i+ is analytic for all rrir, and the 
evolution is just a crossover. 

Next, we discuss the stability of uniform superfluidity 
using two criteria [l^, [H, [2^ [s^, HH : positive definite 
compressibility n(T) and superfluid density p{T) matri- 
ces. 
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FIG. 1: Plots of |Ao|, /i+ and /i_ (in units of eF,+) as a 
function of rrir, when P — 0.5 and l/(fcF,+aF) = (unitarity 
limit). Notice the presence of small cusps in |Ao| and 
when rrir ~ 0.4 which signal a topological quantum phase 
transition discussed below. These cusps are more pronounced 
for higher \P\ (not shown), but are best seen in this figure at 
small grazing angles. 

A. Stability of uniform superfluidity 

In order to analyze the phase diagram at T = 0, we 
solve the saddle point self-consistency (order parameter 
and number) equations for all P and to^ for a set of 
l/(/cF,+aF), and check the stability of saddle point solu- 
tions for the uniform superfluid phase using two criteria. 

The first criterion requires that the compressibility ma- 
trix k{T) is positive definite, where the elements of k{T) 
are Ka,a'{T) = —d'^^a/ [d^iad^ia'). This criterion is re- 
lated (identical) to the condition that the curvature 

of the saddle point thermodynamic potential fJo with 
respect to the saddle point parameter Aq needs to be 
positive Hi, [in. Here, 3^k,± = (3^k,T ± 3^ka)/2 with 
3^k,s = sech^(/3i?k.s/2). Notice that, at low temperatures 
T « 0, 5(i?k,s) — lim^-too /33^k.s/4 where 5{x) is the delta 
function. When at least one of the eigenvalues of k(T), 
or the curvature 9^r2o/5AQ is negative, the uniform sad- 
dle point solution does not correspond to a minimum of 
ilo, and a non-uniform superfluid phase is favored. 

The second criterion requires that superfluid 
density matrix p{T) to be positive definite, where 
the elements of p{T) are defined by pi,j{T) — 
(l//3)Ep{Tr[G^P(p)M]5,,, + fc,fc,Tr[G^P(p)G«P(p)]}, 
where M is a diagonal mass matrix with elements 
Mi.j = ja-nT'o-Sij with 7| = 1 and 7j = —1, and Sij 
is the Kronecker delta. After the evaluation of the 
fermionic Matsubara frequency, we obtain 

p,,iT) - {m^N^ + rrnNi)S,,, - ^ ^ fc.%3^k,+. (17) 

k 
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When at least one of the eigenvalues of p{T) is negative, 
a spontaneously generated gradient of the phase of the 
order parameter appears, leading to a non-uniform su- 
perfluid phase. Notice that the p{T) matrix is reducible 
to the scalar 



Po{T) = m^N^ + m^N^ - ^ ^ k'y^, 

k 



+ 1 



(18) 



in the s-wave case. 

Before discussing ground state phase diagrams, we 
would like to add an additional criterion to fine tune 
the classification of the various phases that emerge as 
a result of unequal masses, interactions, and population 
imbalance. For this purpose we discuss next, the quasi- 
particle excitation spectrum and its connection to topo- 
logical quantum phase transitions. 



B. Topological quantum phase transitions 

The excitation spectrum of quasiparticles is deter- 
mined by energies E'k,! and £'k,2 defined in Eq. Us- 
ing these relations, one can identify surfaces in momen- 
tum space where these energies have zeros, indicating 
that the quasiparticle excitation spectrum changes from 
a gapped to a gapless phase, with a corresponding change 
in the momentum distribution as well. These changes in 
the Fermi surfaces of quasiparticles are topological in na- 
ture. Thus, we identify topological quantum phase tran- 
sitions associated with the disappearance or appearance 
of momentum space regions of zero quasiparticle ener- 
gies when either l/lkp^^ap), P, and/or rrir is changed. 
These topological transitions are shown as dotted lines in 
Figs. [3] through [5l These phases are characterized by the 
number of zeros of i^k.i and ii^k,2 (zero energy surfaces 
in momentum space) such that I) -Ek.i has no zeros and 
i?k,2 has only one, and II) E\^ i has no zeros and £'k,2 has 
two zeros. In the case of unequal masses, the superfluid 
has gapless quasiparticle excitations in both phases I and 
II, when there is population imbalance, e.g., iVf ^ N^. 
We illustrate these cases in Figs, ^a) and[D(b), respec- 
tively, for A^i > Ni. The case of no population imbalance 
P = {Ni = Ni) corresponds to case III, where E'k,! and 
-£'k,2 have no zeros and are always positive, and thus the 
superfluid has always gapped quasiparticle excitations. 

The transitions among phases I, II and III indicate a 
change in topology in the lowest quasiparticle band, sim- 
ilar to the Lifshitz transition in ordinary metals |32l | and 
non-swave superfluids [H, [s^, [H, [H, [13] ■ The topo- 
logical transition here is unique, because it involves an 
s-wave superfluid, and could be potentially observed for 
the first time through the measurement of the momen- 
tum distribution or thermodynamic properties. Notice 
that the topological transition occurs without changing 
the symmetry of the order parameter as the Landau clas- 
sification demands for ordinary phase transitions. How- 
ever, thermodynamic signatures of the topological transi- 
tion are present at low temperatures in the compressibil- 



ity, specific heat, superfluid density, etc., because the 
quasiparticle excitation spectrum changes dramatically. 
The temperature dependence of the quasiparticle contri- 
butions to these properties are exponentially activated 
[~ exp(— _Eg/T)] for the gapped phase (III), or have a 
power law behavior (with different powers of T) in the 
gapless phases (I and II). 
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FIG. 2: Schematic plots of iJk.i (thin lines) and iJk,2 (thick 
lines) versus k for (a) phase I and (b) phase II. Correspond- 
ing momentum distributions rik.T (thin lines) and rik.j, (thick 
lines) are shown in (c) phase I and (d) phase II. Notice that 
the zero of nk,T is slightly shifted upwards from the zero of 
nk.|, for better visualization. 



The zero temperature momentum distributions 

1 - Tff-^k,- Ck. 



2 



2E, 



k.H 



(19) 



for phases I and II are extracted from Eq. ([13]). For mo- 
mentum space regions where ii^k,i > and -Ek,2 > 0, 
the corresponding momentum distributions are equal 
"■k,! = '^k,x- However, when i?k,i > and -Bk,2 < 0, 
then rik.i — and rikj = 1. We illustrate these cases in 
Figs.HKc) and[2Kd) for parameters of Figs, [^a) and[2Kb), 
respectively. Notice that the zero of nk,f is shifted 
slightly upwards to distinguish it from nk.j, in the re- 
gions of momentum space where nk,| = ?^k,j,- Although 
this topological transition is quantum (T = 0) in nature, 
signatures of the transition should still be observed at 
finite temperatures within the quantum critical region, 
where the momentum distributions are smeared out due 
to thermal effects. Although the primary signature of 
this topological transition is seen in the momentum dis- 
tribution, the isentropic ks or isothermal kt compress- 
ibilities and the speed of sound Cg would have a cusp at 
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the topological transition line similar to that encountered 
in |Ao| (see Fig. [ij as a function of the mass anisotropy 
rrir. The cusp (discontinuous change in slope) in K5, kt 
or Cs gets larger with increasing population imbalance. 

Having discussed the finer topological classification 
of possible superfluid phases, we take all the criteria 
together (positive compressibility, positive curvature of 
thermodynamic potential, positive superfluid density, 
and topological character) to present next the resulting 
ground state phase diagrams. 



C. Ground state phase diagrams 

Based on all the previous criteria, we construct the 
P versus phase diagram for seven sets of interaction 
strengths: l/{kF,+aF) = -2, -1 and -0.25 on the BCS 
side shown in Fig. [31 l/{kF^-i.aF) = at unitarity shown 
in Fig. H and l/ikp^+aF) = 0.25, 1 and 2 on the EEC 
side shown in Fig. \5\ In these diagrams, the f (1) label 
always corresponds to lighter (heavier) mass such that 
lighter (heavier) fermions are in excess when P > {P < 
0). Notice that this choice spans all possible population 
imbalances and mass ratios. 

In Figs.IMSl we indicate the regions of normal (N), and 
uniform (U) or non-uniform (NU) superfluid phases. The 
black squares indicate the transition line that separates 
topological phases I and II. In all phase diagrams, phase 
I (II) always appears to the left (right) of the dotted lines 
for P > 0, while phase I (II) always appears to the right 
(left) of the dotted lines for P < 0. 

The normal phase is characterized by a vanishing or- 
der parameter (Aq — 0), while the uniform superfluid 
phase is characterized by po{0) > and d^Qo/dA^ > 0. 
The non-uniform superfluid phase is characterized by 
po{0) < and/or d^flo/dAl < 0, and it should be of the 
LOFF-type havin g on e wavevector modulation only near 
the BCS limit (38l.T39j. although closer to unitarity, we ex- 
pect the non-uniform phase to be substantially different 
from the LOFF phases either having spatial modulation 
that would encompass several wavevectors or presenting 
complete phase separation between paired and unpaired 
fermions. However, from numerical calculations, the su- 
perfluid density criterion seems to be weaker for all pa- 
rameter space and the non-uniform superfluid phase is 
characterized by O^^Iq/OAq < 0, which indicates possible 
phase separation, since at least one of the eigenvalues of 
the compressibility matrix Ko-.o-' also becomes negative. 
Therefore, for a homogeneous system paired (or super- 
fluid) fermions and unpaired (or excess) fermions coexist 
in the uniform superfluid regions, they are phase sepa- 
rated in non-uniform superfluid regions, and the topo- 
logical transition from phase I to phase II may not be 
accessible. However, in a harmonic trap with a large 
superfluid region at the center, the topological phases 
should be observable since the central region is essen- 
tially a "uniform superfluid" with the excess fermions at 
the edge [? ]. The peculiar momentum distribution of 
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FIG. 3: Phase diagram of P = (A^i - N^)/{N^ + N^) versus 
rrir = m-f/mi on the BCS side when a) l/{kF.+aF) = —2 b) 
l/{kF,+aF) = —1 and c) l/{kF,+aF) = —0.25. We show nor- 
mal (N), uniform (U) or non- uniform (NU) superfluid phases. 
The dotted line (black squares) separate topologically distinct 
regions. 



different topological phases would be smeared out by the 
trapping potential, but their marked signatures should 
still be present. Furthermore, these topological phases 
may be accessible in trapped systems at finite tempera- 
ture [4^, or in optical lattices [4l[. 

As shown in Fig.lS^a) and[3I^b), we find a small region 
of uniform superfluidity on the BCS side only when the 
mass anisotropy is small and the lighter fermions are in 
excess {P > 0). Thus, to probe the largest amount of 
phases on the BCS side, mixtures consisting of ^Li and 
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{rrir « 0.15) or ^Li and **^Sr (m^ « 0.07) are good 
candidates. In the rest of the phase diagram, we find a 
quantum phase transition from the non-uniform super- 
fluid to the normal phase beyond a critical population 
imbalance for both positive and negative P. The phase 
space of uniform superfluidity expands while that of the 
normal phase shrinks with increasing interaction strength 
as shown in Figs.[3]Jb) and[2]^c). 
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FIG. 4: Phase diagram of P = (A''t - Ni)/{N^ + Ni) versus 
rrir = /mi at the unitarity limit when l/{kF,+aF) ~ 0. We 
show normal (N), uniform (U) or non-uniform (NU) superfluid 
phases. The dotted line (black squares) separate topologically 
distinct regions. 

This general trend continues into the unitarity limit 
[l/(fci?,+ai?) — 0] as shown in Fig. 21 Since this limit 
is theoreticaUy important as well as experimentally ac- 
cessible, it is useful to analyze the phase diagram as a 
function of population imbalance and mass anisotropy. 
Notice that Fermi mixtures corresponding to mass ranges 
< TO^ < 0.23, like ^Li and ^^Sr {rrir « 0.07) and ^Li and 
^'^K {rrir ~ 0.15) have phase diagrams which are qualita- 
tively different from those corresponding to mass ratios 
0.23 <mr < 0.45 like '^Li and ^^Mg {rrir ~ 0.24) and ^Li 
and {rrir ~ 0.33), since a topological transition line 
may be accessible in the second range. Furthermore, only 
NU and N phases are accessible at unitarity for Fermi 
mixtures in the range of mass ratios 0.45 < rrir < 1 
like and *^Sr {rrir ~ 0.64) or any equal mass mix- 
tures. Notice that our results for the case of equal masses 
{nir = 1) are in close agreement with recent MIT experi- 
ments [J| in a trap. At unitarity, our non-uniform super- 
fluid to normal state boundary occurs at P « ±0.73, and 
the MIT group obtains P w ±0.70(4) for their superfluid 
to normal boundary. 

Additional increase of interaction strength beyond uni- 
tarity on the BEC side leads to further expansion (shrink- 
age) of the uniform superfluid (normal) region as shown 
in Fig. ini^a) and[Sl^b). When heavier fermions are in ex- 
cess {P < 0), a uniform superfluid phase is not possi- 
ble for any mass anisotropy until a critical interaction 
strength is reached. The critical interaction strength 
corresponds to l/{kF,+aF) ~ 0.8 for rrir — 1. Further 
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FIG. 5: Phase diagram of P = (A^i - Ni)/{N^ + N^) versus 
nir = m^/mi on the BEC side when a) 1/ {kp^+ap) — 0.25, 
b) l/{kF,+aF) ~ 1 and c) l/{kF,+aF) = 2. We show nor- 
mal (N), uniform (U) or non- uniform (NU) superfluid phases. 
The dotted line (black squares) separate topologically distinct 
regions. 



increase of interaction strength towards the BEC limit 
[l/(fcF,+aF) > 1], leads to further expansion (shrinkage) 
of the uniform (non-uniform) superfluid region as shown 
in Fig. EI^c), and only the uniform superfluid phase ex- 
ists in the extreme BEC limit [l/{kp^^ap) 3> 1] even for 
P < (not shown). 

Having discussed the ground state phase diagrams, we 
present next fluctuation effects beyond the saddle point 
approximation. 
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IV. GAUSSIAN FLUCTUATIONS 

In this section, we discuss the (gaussian) fluctuation ef- 
fects around the saddle point solutions at finite and zero 
temperatures. Near the critical temperature (T « T^) we 
discuss the time-dependent Ginzburg-Landau (TDGL) 
equation, and at zero temperature (T = 0) we analyse 
the collective phase (or Goldstone) mode as well as the 
effects of harmonic trap in the BEG limit, which are dis- 
cussed next. 



A. Time-dependent Ginzburg-Landau equation 
near critical temperatures 

Our basic motivation here is to investigate the low fre- 
quency and long wavelength behavior of the order pa- 
rameter near Tc where Aq = 0, and derive the TDGL 
equation We use the small q and iv£ ^ uj + i6 ex- 
pansion of 



1 - ?^/(Ck+q/2,T) - "/(Ck-^q/2,i) 

k 



FkP, 



Ck+q/2,T + Sk-q/2,i ^ 

where L~^{q) ~ F^\{q), to obtain the TDGL equation 



A{x) = (20) 



in the real space x — (x, t) representation. Here, nf[x) — 
l/[exp{(jx) + 1] is the Fermi distribution. 

Expressions for the coefhcients a, 6, Cij and d are pre- 
sented in appendix [Bl The condition a = corresponds 
to the Thouless criterion, and the coefficient of the non- 
linear term is positive (6 > 0) guaranteeing the stability 
of the effective theory. The kinetic energy coefficient is 
an effective inverse mass tensor which reduces to a scalar 
c in the s-wave case. The time-dependent coefficient d 
is a complex number, and its imaginary part reflects the 
decay of Cooper pairs into the two-particle continuum for 
> 0. However, for /x+ < 0, the imaginary part of d 
vanishes and the behavior of the order parameter A(a;) is 
propagating reflecting the presence of stable (long lived) 
bound states. 

Since a uniform superfluid phase is more stable in the 
BEG side, we calculate analytically all coefficients in 
the BEG limit where |/x±| |ef,|/2 > Tc. We obtain 
a = oi + 02 = —Vm\{2^^ — eb)aF /{8tt) + Vm+nea'^, 
b ^ bl + b2 ^ Fm3.a|,/(167r) - Vml{dne/dfie)a%, 
c — Vm'^ap /[8Tr{m-i -|- mj^)], and d — Vrn^ap/iS''!^)- 
Here, e labels the excess type of fcrmions and Ue is 
the density of unpaired fermions. Through the rescaling 
^{x) = ^/dA{x), we obtain the equation of motion for a 
dilute mixture of weakly interacting bosons and fermions 



fiB'i'ix) + [Ubb\'^{x)\'^ + UBFne{x)]'i'{x) 



\/^^{x) .9*(j;) 



2m B 



dt 



= 0, 



(21) 



with bosonic chemical potential (Ib = —ai/d = 2^+ — Cb, 
mass iriB — d/c = + m^, and repulsive boson- 
boson Ubb — Vbi/d^ ~ ATTap/rn^ and boson-fermion 
Ubf = Vai/{dne) = 8TraF/'m+ interactions. This pro- 
cedure also yields the spatial density of unpaired fermions 
given by 



neix) = [a2/d + b2\-^ix)\yd^]/UBF 

= ne-[/BF(ane/a/ie)|*(x)p. 



(22) 



Since dng^/dfj,g > the unpaired fermions avoid regions 
where the boson field |^'(a;)| is large. Thus, in a harmonic 
trap, the bosons condense at the center and the unpaired 
fermions tend to be at the edges. Notice that, Eq. ([2T|) 
reduces to the Gross-Pitaevskii equation for equal masses 
with P = 3, and to the equation of motion for equal 
masses with P ^ [13. 

Next, we recall the standard definitions of the in- 
teractions in terms of the scattering lenghts Ubb — 
Attubb /"rnB and Ubf ^ Attubf /mBF, where tob is the 
mass of the bosons and uibf — '^memB /{rriB + We) 
is twice the reduced mass of a boson of mass tub and 
an excess fermion of mass rUe- Gombining these defini- 
tions with our results for Ub b and Ubf va. terms of the 
fermion- fermion scattering length of, we can directly re- 
late the boson-boson scattering parameter 



O^BB 



niB 
1 



-ap 



1 



mi 



2m 



2m 



ap 



and the boson-fermion scattering parameter 



qbf 



-ap 



(23) 



(24) 



m+{mB + 

to ap. Notice that these expressions reduce to ubb = 
2aF and aBP — 8ap/3 for equal masses [1, [l3l- A better 
estimate for aBB can be found in the literature [43| . 

Since the effective boson-fermion system is 
weakly interacting, the BEG temperature is 
Tc = 27r[r7.s/C(3/2)]2/3/ms, where C{x) is the Zeta 
function and ub = (n ~ ne)/2. 

Notice that, the effective total number equation for the 
boson-fermion mixture can be written as 



5I"/(Ck,e)+2 



q 



rib 



2mB 



fJ'B 



(25) 



where rib(e) = l/[exp(/3e) — 1] is the Bose distribution 
and JiB — > 0^ includes also the Hartree shift. In the limit 
when Tc — > 0, we obtain the critical chemical potential for 
unpaired fermions at the normal-to-stable uniform super- 
fluid boundary as given by /Xe = 2^/'^(m+/me)eF,-i-, where 
e = (T7 i) labels excess type of atoms. Since, /Lt+ ef,/2 
in this limit, the critical chemical potential imbalance 
is given by /i_ = 7e[— eb/2 -I- 2^/'^(m+/me)eF,+], where 
tb = — l/(m+a|,) is the binding energy, and 71 = +1 and 

7i = -1- 

This concludes our analysis for the homogenous mix- 
ture of two types of ultra-cold fermions at flnite tem- 
peratures. Next we discuss collective excitations at zero 
temperature. 
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B. Sound velocity at zero temperature 

In order to obtain the collective mode spectrum, we 
use the effective action defined in Eq. ^ and express 
A{q) = [\{q) + i9{q)]/V^ in terms of the amplitude \{q) 
and phase 9{q) fields, respectively. Using the matrix ele- 
ments of defined in Eqs. ([5]) and © and described 
in appendix [XI we can obtain the matrix elements of the 
fluctuation matrix in the rotated basis {X{q),6{q)). The 
diagonal elements of the fluctuation matrix in the rotated 



basis become M.^\{q) — [F^ } - 
and M,-j(g) = [Frj - F^^ - F 
off-diagonal elements become M 



F2',2]/2 

+ F^lil/i, while the 



1,2 



F^']/2 with M^\{q) = (M7i)*(g) 



The collective modes are found from the poles of the 
fluctuation matrix 'M.(q) determined by the condition 
detM^^((7) = 0, when the usual analytic continuation 
iv£ w + iO~^ is performed. The easiest way to get the 
phase collective modes is to integrate out the amplitude 
flelds to obtain a phase-only effective action. To obtain 
the long wavelength dispersions for the collective modes 
at r = 0, we consider |P| ^ or kF,+ — kp,^ = kF,i 
limit, and expand the matrix elements of F~^{q) to sec- 
ond order in |q| and w to get 'M'^\{q) = A+C\q\'^ — Dw^ , 

^eeil) = Ql^l^ ^ ^'^^ M^g(g) = iBw, such that 



A- 



—iBw 



iBw 
Q\q\^-Rw^ 



The expansion coefficients are given in the Appendix |B] 
Thus, there are two branches for the collective excita- 
tions, but we focus only on the lowest energy one corre- 
ponding to the Goldstone mode with dispersion w(q) = 
w|q|, where 



AQ 



AR + B"^ 



(26) 



is the speed of sound. Extra care is required when P ^ 
since Landau damping causes collective excitations to de- 
cay into the two-quasiparticle continuum even for the s- 
wave case, since gapless fermionic (quasiparticle) excita- 
tions are present (see Fig. [2]). 

The BCS limit is characterized by the criteria /i+ > 
and /i+ w eF.+ 3> |Ao|. The expansion of the ma- 
trix elements to order |qp and w'^ is performed under 
the condition [w, |qp/(2TO+)] <C |Ao|. The coefficient 
that couples phase and amplitude flelds vanish {B — 0) 
in this limit. Thus, there is no mixing between the 
phase and amplitude modes. The zeroth order coeffi- 
cient is A ^ V, and the second order coefficients are 
C = 0/3 = VvF,'\VF,i/ i5G\Ao\'^), and D = i?/3 = 
X>/(12|Aop). Here, VF,a = kF,a/m„ is the Fermi velocity 
and V = m+VkF,+/{'2n'^) is the density of states per spin 
at the Fermi energy. Thus, we obtain 



(27) 



with Vcr = VF,a/V3, which reduces to the Anderson- 
Bogoliubov relation when the masses are equal. 



0.6 



0.4 



0.2 




0.25 



0.5 



0.75 



FIG. 6: Sound velocity v (in units of vf.+ ~ kF, + /m+) ver- 
sus rrir for \/{kF,+aF) = —1 (solid line), l/(fcF,+csF) = 
(solid squares) and l/{kF,+aF) ~ 1 (hollow squares). Here, 
populations are equal (P = 0). 

On the other hand, the BEG limit is characterized by 
the criteria /i+ < and ^k,+ 3> |Ao|. The expansion of 
the matrix elements to order |qp and w'^ is performed 
under the condition [w, |qp/(2m+)] <^ |/^_|_|. The coef- 
ficient B indicates that the amplitude and phase 
fields are mixed. The zeroth order coefficient is ^ = 
k|AoP/(2|/j,_|_|), the first order coefficient is B = k, and 
the second order coefficients are C = Q = K/[2{mi+mi)] 
and Z) = i? = k/(8|;U+|), where k = V/i32^\fi+\eF,+). 
Thus, we obtain 



lAn 



^4(m| +mj|/i+| 



(28) 



with Va- = ^ ^TiUiyaF I iv^c ■ Notice that the sound velocity 
is very small and its smallness is controlled by the scat- 
tering length aF- Furthermore, in the theory of weakly 
interacting dilute Bose ga s, the sound velocity is given 
by = ^ A.'KaBB'n.B I '^B- Making the identification that 
the density of pairs is ns — n+ , the mass of the bound 
pairs is tob = + mi and that the Bose scattering 
length is 



O-BB 



niB 



-gf = 



1 



mi 



2m 



2mx 



aF, 



(29) 



vb reduces to the well known Bogoliubov relation when 
the masses are equal. Therefore, the strongly interact- 
ing Fermi gas with two species can be described as a 
weakly interacting Bose gas at zero temperature as well 
as at finite temperatures [l^]. Notice that ubb reduces 
to aBB — 2aF for equal masses [l^l in the Born approx- 
imation. A better estimate for ubb can be found in the 
literature j43| . 

In Fig. [51 we show the sound velocity as a function of 
the mass ratio for three values of the scattering pa- 
rameter l/{kF,+aF) = —1)0 and 1 corresponding to the 



9 



BCS side [l/{kp^+aF) — —1], unitarity [l/^kp^+ap) ~ 0], 
and to the BEC side [1/ (kp^^ap) — !]■ Notice that 
the speed of sound could be measured for a given rrir 
using similar techniques as in the single species case 
rrir = 1 Hlil. 

This concludes our analysis for collective excitations 
in the homogenous mixture of a two types of fermions 
at zero temperature. We discuss next the effective field 
theory between paired and unpaired fermions. 



C. Weakly interacting paired and excess fermions 
(Bose-Fermi mixtures) at zero temperature 

In this section, we concentrate on the BEC regime, 
where the paired and unpaired (excess) fermions can be 
described by a mixture of molecular bosons and fermions. 
In this limit, the resulting equation of motion is identical 
to Eq. (pij) near the critical temperature, except that 
all parameters are evaluated at zero temperature. Thus, 
at low temperatures the system continues to behave as 
a dilute mixture of weakly interacting bosons (formed 
from paired fermions) and unpaired fermions, and can 
be described by the free energy density 

.F(x) = E{^) ~ fi,n,{^) - Aii3|*(x)|^ (30) 

where the energy density is 

^(x) =Kb+Kp + ic/BB|^'(x)|4 + C/BFrJe(x)|*(x)|2. 

Here, Kb is the kinetic energy density of bosons (as- 
sumed to be much smaller than all the other energies) 
and Kp is the kinetic energy density of fermions. Aver- 
aging these energy densities over the spatial coordinates 
F — J dx!F{x)/V leads to a ground state average free 
energy density 

1 , 3 

F = -UbbTIb + UBpneUB + -ep.ene - fJ-ene - ^J-snB, 

2 5 

where Ue {ub) is the average density of fermions (bosons), 
and ep,e is the fermi energy of the excess fermions. Using 
the positivity of the Bose-Fermi compressibility matrix 
Ka,f3 = d^a/dufj, whcrc a,/3 = e, -B, one can show that 
bosons and fermions phase separate when the condition 



Up > 



47r 



TT 

nie 



(Ubb\ 
Wbf) 



(31) 



is satisfied (46|. 

Using the boson-boson and boson-fermion interactions 
Ubb = '^T^a.BB /rrLB and Ubf ^ 4:TraBP /msp, the scat- 
tering paramaters indicated in Eqs. and ([M]) . as well 
as the relations |P| = rie/n, n = n-^ + = fc|,^/(37r^) 
and riB = (n — ne)/2, then phase separation occurs when 



IPI > 



-f-V 



kp^j^ap 



(32) 



Here e labels excess type of fermions, and m,, is the 
mass of the unpaired fermions, and rrij^. is twice the 
reduced mass of the t and | fermions. This expres- 
sion is quantitatively correct in its region of validity, 
i.e., when l/^kp^+ap) ^ 1, however, it still gives semi- 
quantitative results for l/{kp^^ap) 2. For instance, 
in the case of an equal mass mixture, this expression 
would suggest that the resulting Bose-Fermi mixture is 
uniform when l/{kp_^aF) > 2.5 for |P| — *■ 0.5, and when 
l/{kp^^ap) > 3.2 for \P\ 1. From our numerical cal- 
culations we find l/{kp^^ap) > 1.9 for \P\ 0.5, and 
when l/{kp^+ap) > 2.4 for \P\ 1. 

The analytic expression given in Eq. (|32p may be used 
as a guide for the boundary between phase separation 
(non-uniform) and the mixed phase (uniform) for any 
mixture of fermions. In particular, this relation serves 
as an estimator for the phase bounday for future ex- 
periments performed in the BEC limit of unequal mass 
fermions with population imbalance. This relation can 
be rewritten in terms of the mass ratio rrir = fn]' /^j, by 
realizing that the ratio m+/me = 2to,./(1 -|- rrir) when 
I (heavier) fermions are in excess, and that m^/rrie — 
2/(1 -|-mr) when "f (lighter) fermions are in excess. Thus, 
when I (heavier) fermions are in excess, the critical po- 
larization below which phase separation (non-uniform 
phase) occurs is 



p(i) _ 

^c,l - - 



1 



2 \4:kp^+ap 



1 



(33) 



while, when | (lighter) fermions are in excess, the critical 
polarization above which phase separation (non-uniform 
phase) occurs is 



2 \4:kp,-^-ap 



1 



(34) 



(1) 



Notice that in the equal mass {rrir — 1) case P, 

— 1^ as required by symmetry. 

In addition, we can also describe analytically a finer 
structure of non-uniform (phase separated) superfluid 
phases deep into the BEC regime. For a weakly inter- 
acting Bose-Fermi mixture, the phase separated region 
consists two regions: PS(1), where there is phase sep- 
aration between pure fermions and pure bosons (tightly 
paired fermions), and PS (2), where there is phase separa- 
tion between pure fermions and a mixture of fermions and 
boso ns ( tighly bound fermions). Following the method of 
Ref. [401, we obtain analytically the condition 



rip > 



11257r4 U 



BB 



128m4 



:nB, 



(35) 



for the transition from the PS(2) to the PS(1) phase. 

Using our effective boson-boson (Ubb) and effective 
boson-fermion Ubf interactions, we can rewrite this re- 
lation as 



\P\ > — 
' ' - 11 



157r 



/rrie 
kp^j^ap 



- ^, (36) 
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FIG. 7: Phase diagram of \P\ = |7V| - Ni\/{Ni + iV|) versus 
{m+frrie) /{kF,+aF) in the BEG Hmit. We show the uniform 
superfluid (U) phase where paired and unpaired fermions co- 
exist, and the phase separated non-uniform superfluid phases 
PS(1) and PS(2). The PS(1) region labels phase separa- 
tion between pure unpaired (excess) and pure tightly paired 
fermions (bosons), while the PS (2) region labels phase sepa- 
ration between pure unpaired (excess) fermions and a mixture 
of unpaired and tightly paired fermions. 



where we used ns = (rt — ne)/2 as the boson density 
and |P| = rie/n. These phase boundaries can also be ex- 
pressed in terms of rrir and l/{kF,+aF) as follows. When 
I (heavier) fermions are in excess, the critical polarization 
below which the transition from PS(2) to PS(1) occurs is 



p(2) _ _ 1 



157r 



11 V IB/cF. + flF 



1 



11' 



while, when | (lighter) fermions are in excess, the critical 
polarization above which the transition from PS (2) to 
PS(1) occurs is 



d(2) 



1 / IStt 



11 V 16fci? j^ap 



1 



IT- 



Notice that in the equal mass 



1) case P, 



-P' 



(2) 



(2) 



J, I as required by symmetry. 
In Fig. [71 we show phase diagram of uniform and non- 
uniform superfluidity as a function of population imbal- 
ance |P| and (m^ I rrif,) I (kp^j^aF)^ which is strictly valid 
in the BEC limit when Xj (kp j^ap') 3> 1. In this figure, 
we show the uniform superfluid (U) phase where tightly 
paired and unpaired fermions coexist, and phase sepa- 
rated (non-uniform) superfluid (PS) phases. The PS(1) 
region labels phase separation between pure unpaired 
(excess) and pure tightly paired fermions (bosons), while 
the PS (2) region labels phase separation between pure 
unpaired (excess) fermions and a mixture of unpaired and 
tightly paired fermions. The phase boundary between U 
and PS (2) phases is determined from Eq. (|32p. and the 
phase boundary between PS(2) and PS(1) phases is de- 
termined from Eq. (j36p . For a fixed mass anisotropy m^, 
when |P| is large, we find phase transitions from PS(1) to 



PS(2) to U phase as the interaction strength Xj (kp^j^ap) 
increases. However, when |P| is small, we find a phase 
transition from the PS(1) to the U phase as Xjikp^j^ap) 
increases. This last result is unphysical and signals the 
breakdown of the approximations. 
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FIG. 



8: Phase diagram of P = (iVj - iVi)/(iVt + ^i) versus 
rrir — m-^/my in the BEG limit when (a) l/{kF,+aF) — 3 
and (b) l/(A:F,+fflF) ~ 4. We show the uniform superfluid (U) 
phase where paired and unpaired fermion coexist, and phase 
separated non-uniform superfluid phases PS(1) and PS (2). 

In Fig. [U we show the phase diagram of popula- 
tion imbalance P versus mass anisotropy = TO| /mi 
in the BEC limit when (a) l/{kp^+ap) = 3 and (b) 
l/{kp^^ap) — 4. We indicate the uniform superfluid (U) 
phase where paired and unpaired fermion coexist, and 
phase separated non-uniform superfluid phases PS(1) 
and PS(2). The phase boundary between U and PS(2) 
phases is determined from Eq. p3p when P < 0, and 
from Eq. ([M)) when P > 0. In addition, the phase bound- 
ary between PS (2) and PS(1) phases is determined from 
Eq. ([27]) when P < 0, and PS(1) phase does not ex- 
ist when P > 0. Notice that these phase diagrams are 
very similar to the one given in Fig. [HJc) , with the added 
refinement of the non-uniform superfluid phases PS(1) 
and PS(2). For a fixed interaction strength 1/ {kp,-^-ap), 
when |P| is large, we find phase transitions from PS(1) 
to PS (2) to U phase as the mass anisotropy nir increases. 
However, when |P| is small, we find a phase transition 
from the PS(1) to the U phase as nir increases. This last 
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result is unphysical and again signals the breakdown of 
the approximations. 

To summarize, we analyzed analytically the structure 
of non-uniform (phase separated) superfluid phases in 
the BEC regime. However to understand the ultracold 
atomic experiments, one needs also to consider the trap- 
ping potential, which is discussed next. 



D. Effects of trapping potential 

For simplicity, we approximate the trapping potential 
by an isotropic harmonic function where the potential 
energy is Vcr(x) — Q;cr|xp/2 such that the local chemical 
potentials are given by 



/io-(x) = l^la - ^Qfcrlxp. 



(39) 



Here, a^r is proportional to the trapping frequency of 
the a type fermion, which is typically different for each 
kind of atom. In general, it is quite difficult to make 
completely isotropic traps and harmonic traps are typi- 
cally elongated such that VJj(x) — a^^xx'^ /2 + aa,yy'^ + 
Cia^zZ^ with cko-^x — Oia.y ^ O-a^z- Howcvcr, the same 
qualitative behavior occurs in the elongated or spher- 
ically symmetric (isotropic) traps, and we will confine 
ourselves for simplicity to the isotropic case. When ex- 
perimental data becomes available and all the numbers 
are known, one can revisit this problem for detailed com- 
parison between theory and experiment. 

Again, we confine our discussion to the BEC regime, 
and obtain the equation of motion for a dilute mixture 
of weakly interacting bosons and fermions at zero tem- 
perature 



V-t(x) + nW]^(^) - ^77^ = (40) 



v2*(.t) .a*(x) 

= t 

2m B dt ' 

where the spatial density of unpaired fermions is 



ne{x) 



= (x) 



Ub.'-^Mx)?. 



(41) 



These results are quite similar to the case of equal 
masses [l3l. Notice that setting Va = Q reduces the 
problem to the free space case discussed in the previous 
subsection. Here, ne(x) — (1/V^) X^k [^k,e — /^e(x)] , 
where n/(e) — l/[exp(/3e) -I- 1] is the Fermi distribution. 
In the BEC limit when ap — > 0^, we can approximate 
the local density of unpaired fermions as 



ne{x) 



[ek,e-/ie(x) + C/si.|*(x)|2] , (42) 



which at zero temperature leads to 



1 



ne(x) = g^{2me [/Xe(x) - UBFUBi^W^- (43) 
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FIG. 9: (a) Density Tio-dxl) of fermions (in units of k%^^), and 
(b) total density of fermions n(|x|) (hollow circles) and density 
of unpaired fermions ne(|x|) (solid circles) versus trap radius 
|x|/|x|tf. Here P = 0.5 and l/{kF,+aF) = 2. In (b), we also 
compare w(|x|) for P = when l/{kF,+aF) = 2 (crosses). 



Notice that the density of bosons at zero temperature 
is given by n_B(x) = |5'(x)p. Therefore, we need to 
solve Eq. (|40l) self-consistently with the number of un- 
paired (excess) N^. ~ J dxnf, (x) and paired (bound) 
N}yf — 2 J rfxnB(x) fermions such that the total num- 
ber of fermions is N = + N^f. 

Next we solve the self-consistency equations for a ^Li 
and mixture (m^ = 0.15) within the Thomas-Fermi 
(TF) approximation, where the kinetic energy term in 
Eq. is neglected. This leads to a coupled equation 
for density of paired and unpaired fermions 



"-b(x) 



fiB - V't(x) - Vi(x) - t/gf ne(x) 
Ubb 



(44) 



In our numerical analysis, we choose for convenience 
ai = ai = a, \/ ikp^-^-ap) — 2 and P ~ 0.5 such 
that — iNi. However, in a more realistic case 
a I ^ a I — a, since atoms with different masses may 
experience different trapping potentials due to their dif- 
ferent polarizabilities. 

In Fig. [IKa), we show the density n^dxl) of a type 
fermion (in units of A:|,^_) as a function of |x|/|x|tf, 
where |x|tf is the TF radius defined by cf.+ = 
kp ^/{'2m+) — a\x\'^p/2. We also scale the total num- 



12 



ber of fermions with TV = fc|,_|.|x||,^/24. We find that 
the density of "f and | fermions are similar close to the 
center of the trapping potential, while most of the excess 
fermions are close to the edges. In Fig. Eljb), we show 
the total density n(|x|) = ?T.-f(|x|) + nx(|x|) as well as 
unpaired fermions ne(|x|) = n|(|x|) — n|(|x|). In both 
figures, we find a clear indication of phase separation be- 
tween paired and unpaired fermions. In Fig.lHfb), we also 
compare the total density of fermions JT-dxl) for the same 
parameters when the populations are balanced Nf = . 
When P 7^ 0, the total local density of fermions at the 
center of the trap is reduced in comparison to the P = 
case for the same fermion scattering parameter, since 
the unpaired fermions are pushed away from the cen- 
ter of the trap due to Ubf- These findings for unequal 
masses are similar to previous results on equal mass mix- 
tures [liHIiliiiil,!!!. 

Lastly, for the parameters used, the harmonic trap 
tends to favor phase separation into a PS(l)-type phase 
where one has almost pure fermion and almost pure bo- 
son regions. In a harmonic trap it may be still possible 
to realize the PS(2)-type phase where one has an almost 
pure fermion region and an almost pure mixed phase of 
bosons and fermions, provided that good control over the 
trapping potentials is possible. Having concluded our dis- 
cussion of the effects of a trapping potential, we present 
next the summary of our conclusions. 



V. CONCLUSIONS 

In summary, we analyzed the phase diagram of ultra 
cold mixtures of two types of fermions (e.g., ^Li and '^'^K; 
^Li and s^Sr; or and ^^Sr) from the BCS to the 
BEC limit as a function of scattering parameter, popula- 
tion imbalance, and mass anisotropy. We found that the 
zero temperature phase diagram of population imbalance 
versus scattering parameter is asymmetric for unequal 
masses, having a larger stability region for uniform su- 
perfluidity when the lighter fermions are in excess. This 
result is in sharp contrast with the symmetric phase di- 
agram for equal masses. 

In addition, we discussed topological quantum phase 
transitions associated with the disappearance or appear- 
ance of momentum space regions of zero quasiparticle en- 
ergies when either the scattering parameter or population 
imbalance are changed. These quantum phase transitions 
are reflected in the momentum distribution as well as in 
thermodynamic properties, however they seem to lie in 
the non-uniform region of the phase diagram, but may 
survive at the center of a harmonic trap (37j . Further- 
more this phase may be observable at finite temperatures 
in trapped systems [i^, or in optical lattices 41]. 

We also analyzed gaussian fluctuations around the sad- 
dle point order parameter both at finite and zero tem- 
peratures. Near the critical temperature, we derived 
the Ginzburg-Landau equation, and showed that it de- 
scribes a dilute mixture of composite bosons (tightly 



bound fermions) and excess (unpaired) fermions in the 
BEC limit. At zero temperature, we obtained analyti- 
cally the dispersion of collective excitations in the BCS 
and BEC limits, and showed numerically the evolution 
from the BCS to BEC regimes in the case of zero popu- 
lation imbalance. In addition, we discussed analytically 
how phase separation between paired fermions and excess 
fermions emerges analytically at zero temperature in the 
BEC limit. Lastly, we discussed the effects of a harmonic 
trapping potential, and concluded that phase separation 
between paired and unpaired fermions is favored even in 
the BEC limit. 
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APPENDIX A: INVERSE FLUCTUATION 
PROPAGATOR 

In this Appendix, we present explicitly the elements of 
the inverse fiuctuation propagator F~^{q). The diagonal 
matrix element of F^^(q) is given by 

i_l 2 ^fiE^,-Ki)-ME.^Ki) 
-•^11^" + / 1 ^'a-k^a+k — : \ — -p, -p, 



'a-k'^a+k 



,2 2 



n/(^q_k,2) " '^/(-E2+k,2) 

ivi - i?^_k,2 + -Sa+k,2 
1 - ri/(£;a+k4) - "/(-^t-k,2) 



}|r: 



2 1 - "/(^f-k.l) - '^/(-Ba+k,2) 
k«a_k- 

|2 



(Al) 



and the off-diagonal matrix element of F ^{q) is given 
by 



^,2 = X]"f+k"f+k^f-k^f-k{ 
k 

"/(-^a-kg) -"/(-g2+k,i) 
ivi + Pa-k,! - £'f +k,i 

n/(£;|_k,2) -"/(-ga+k,2) 
ivi - £'a_k,2 + £^a+k,2 

1 - ri/(Sa+k,i) - n/(£;q_k^2) 



k.2 



+ 



if^ + Sq_k,l+£;q+k,2 
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(1 + a,+/i^k,+)/2 and 



(1 - 



l/[exp(/?x) + 1] is the Fermi 



where 

Ck,+/£^k,+)/2, and 
distribution. 

For the s-wave case considered in this manuscript, a 
well defined low frequency and long wavelength expan- 
sion is possible in two limits: (I) at zero temperature 
(T = 0) when population imbalance is zero P = such 
that the Fermi functions in Eqs. (jAip and (|A2p vanish, 
and (11) near the critical temperature (T « T^) where 
|Ao| — > such that waj-k (when ^s±k + > 0) "3±k 
(when Ca±k.+ < 0) in Eqs. (IXTj) and (IX2|) vanish. Other 
than these two limits, there is Landau damping which 
causes the collective modes to decay into the two quasi- 
particle continuum. 



APPENDIX C: EXPANSION COEFFICIENTS AT 

T = 

In this Appendix, we perform a small q and ivi —f 
w + iO~^ expansion of the effective action at zero temper- 
ature {T — Q). From the rotated fluctuation matrix 
expressed in the amplitude-phase basis, we can obtain 
the expansion coefficients necessary to calculate the col- 
lective modes. We calculate the coefficients only for the 
case of zero population imbalance P = 0, as extra care is 
needed when P 7^ due to Landau damping. In the long 
wavelength (|q| 0), and low frequency {w 0) limits 
the condition {w, <C min{2£'k,+}, is used. 

The coefficients necessary to obtain the matrix element 
M^\(g) are 



APPENDIX B: EXPANSION COEFFICIENTS AT 

T = Tc 

In this Appendix, we derive the coefficients a,b,Cij,d 
of the time dependent Ginzburg-Landau theory described 
in Eq. ((20|) . We perform a small q and iv^ w + iO+ 
expansion of the effective action near the critical tem- 
perature (T « Tc), where we assume that the fluctuation 
field A(x, t) is a slowly varying function of x and t. The 
zeroth order coefficient L~^(0, 0) is given by 



E 



-^k,+ 

26^ 



(Bl) 



where Xk,± = (-''^k,! i -'^k,i)/2 and Xk,<T = 
tanh(/9^k,(T/2). The second order coefficient 
d'^L~^{c[,0)/{dqidqj) evaluated at q is given 

by 



E 

k 



k,T^k,T 



Xk iFi 



rn' 



I 



32^^+ 



k^i k j 



2m+ mi^k,+ 



i6ek.+ 

-''^k,+ 



(B2) 



where C± = (lk,T/"^T ^ ^k,i/'^i)/2, and yk,cr — 
sech^(/3^k,cr/2). Here, Sij is the Kronecker delta. The 
coefficient of the fourth order term is approximated by 
its value at q = 0, and given by 



b = 



k,+ 



k 



4^k,+ "'.k,+ 



iFkl 



(B3) 



The time-dependent coefficient has real and imaginary 
parts, and for the s-wave case is given by 



lim V X] 



k,+ 



1 



2w 



irkP 



(B4) 



where 6{x) is the Delta function. 



^-E 



|Ao| 



k 



(CI) 



corresponding to the (q = 0, w = 0) term, 



C 



E 



10|Ao 



^^k,+ -3|Aop 



'^kH 



lOlAol 



|Ao| 



m_ 



(C2) 



corresponding to the |qp term, and 



^ = E 



8^k,+ 



(C3) 



corresponding to the uP' term. 

The coefficients necessary to obtain the matrix element 
M.gUq) are 



Q-E 



6 



^k,+ -3|Aop 



Eh 



8™+^k,+ 



corresponding to the |qp term, and 

p = y^ 
k 



(C4) 
(C5) 



corresponding to the term. 

The coefficient necessary to obtain the matrix element 



6 



k 



(C6) 



corresponding to the w term. These coefficients can be 
evaluated in the BCS and BEC limits, and are given in 
Sec. UVBl 
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